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1.   SUMMARY  AND  INTRODUCTION 

In  IMM-NYU  Report1  No.  226  (195'6) ,  the  equations  governing 
geostrophic  vortex  motion  were  derived  and  the  possibility  of 
applying  the  equations  to  the  problem  of  hurricane  tracking  was 
discussed.  Since  that  time  E.  Isaacson,  D.  A.  Levine  and  the 
author  have  had  some  success  in  the  short  range  (one  or  two 
days)  prediction  of  hurricane  motion  over  the  Atlantic  Ocean  ap- 
proaching the  eastern  seaboard.   In  this  report  we  present  an 
outline  of  the  work  and,  briefly,  what  is  planned  for  the  imme- 
diate future. 

The  basic  equations  of  geostrophic  vortex  motion  derived  in 
GVM  are 

■(A-^2H  ■  0  (la) 


dX, 
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(   )  =  (   )  +u(1).(   )  +v(1).(   )  (lb) 


dt  t 


U(D  =  .^  ,   v(1)  =  1rx  •  do) 

The  third  order  non-linear  differential  equation  (la)  for  the 
stream  function  \|/  =£  gh'  ' /f   is  called  the  Geostrophic  Conserva- 
tion Equation.   Subscripts  denote  partial  differentiation; 
/\(   )=()+()    is  the  two  dimensional  Laplacian  operator; 


Based  on  this  report,  a  manuscript  titled  "Geostrophic  Vortex 
Motion"  has  been  submitted  to  the  Journal  of  Fluid  Mechanics, 
Cambridge  University,  England,  for  publication;  in  this  re- 
port we  refer  to  this  paper  as  GVM. 
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KT  =  f  /gh  where  f (=  constant)  is  the  Coriolis  parameter,  g  is 

the  gravitational  acceleration  and  h  is  the  depth  of  the  homo- 

o 

geneous  atmosphere  at  rest;  and   =  at  is  the  scaled  time  where 
a  «  1  is  a  small  parameter  related  to  the  strength  of  distur- 
bances. The  superscript  (1)  denotes  that  the  solution 
(h*1',^1',^1')  is  the  first  order  approximation  to  the  exact 
solution  (h,u,v)  =  (h0+ah(1^+. . .,  au(1^+...,  av(1'+...)  expres- 
sed as  a  power  series  expansion  with  respect  to  the  parameter  a, 
perturbed  on  the  flow  at  rest  (h,u,v)  =  (hQ,0,0). 

The  singular  vortex  solution  of  (1)  which  vanishes  at  infi- 
nity is 

Vl^o')  =  "  HKo(K-lr-rol>  •  (2a) 

The  vortex  strength  y  is  defined  by 

t  r 

p  M 
Y  =Cj>\|/  ds-  f^.      j  \|/dxdy  (2b) 

C  R 

where  the  integration  is  over  a  region  R  with  boundary  C  which 
encloses  the  vortex  and  v  denotes  differentiation  in  the  outward 
normal  direction  to  G.  K  is  the  modified  Bessel  function  of 


the  second  kind  and  zero  order  of  argument  lt|r-r0|  where 

|r-rj2  =  [x-xo(T)]2+[y-yo(r')]2  (2c) 


(2a)  is  called  the  Geostrophic  Point  Vortex. 

For  application  of  equations  (1)  and  (2)  to  the  problem  of 
hurricane  tracking  we  write  the  stream  function  \|i  as  the  sum  of 
the  singular  solution  (2a),  representing  the  hurricane,  and  the 
remaining  continuous  regular  solution  ^  (see  GVM) : 
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*  =  Vfr-ro,)'H|'l(x'y'r)  •  (3a) 

Putting  (3a)  in  (1)  we  obtain  the  equations  of  motion  of  a 
single  geostrophic  point  vortex  embedded  in  a  continuous  flow 
field  (see  GVM) : 

where 


(3b) 


and 


u(1)  =  -(V^y  '  v(1)  =  ( Wx  (3c) 

d^ 

where  -(^, )_  and  (^-,)__  are  the  x-  and  y-components  of  the  velo- 
X  y        XX 

city  of  the  continuous  flow  field  \|/,  evaluated  at  the  vortex 

Point  U0,y0). 

When  the  initial  and  boundary  conditions  are  given,  either 

the  system  of  equations  (1)  or  (3)  can  be  applied  to  the  problem 

p 

of  hurricane  tracking.   There  is  some  feeling  among  meteorolo- 
gists that,  in  numerical  prediction  of  hurricane  motion,  the 
hurricane  must  be  explicitly  distinguished  from  the  rest  of  the 
flow  as  shown  in  (3a)»  Otherwise  the  hurricane  which  is  a 
strongly  concentrated  disturbance  might  not  be  resolvable  by  the 
usual  spacial  grid  spacing  (one  to  two  hundred  miles).   Our  in- 
vestigations so  far  indicate  that  equations  (1)  can  be  directly 
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In  this  report  we  call  the  application  of  equations  (1)  the 
Total  Flow  Method,  and  the  application  of  equations  (3)  the 
Point  Vortex  Method. 


1+ 

applied  without  losing  the  hurricane  In  the  grid  mesh.  Thus, 
our  present  position  on  this  question  is  that  advantages  and  dis- 
advantages exist  for  either  method  of  hurricane  tracking  and  an 
important  part  of  our  overall  program  is  to  determine  the  condi- 
tions under  which  one  method  is,  if  at  all,  better  than  the 
other.   It  is  rather  apparent  that  using  equations  (1)  entails 
less  numerical  work.  Also,  the  determination  of  preference  of 
one  method  over  the  other  appears  to  be  strongly  dependent  on 
the  reliability  and  accuracy  of  the  initial  data  particularly 
in  the  neighborhood  of  the  eye  of  the  hurricane,  i.e.,  of  the 
order  of  one  hundred  miles  radius.  Thus  we  expect  some  shift  of 
the  relative  advantages  of  one  method  over  the  other  as  better 
initial  data  is  obtained. 

The  remainder  of  this  report  is  divided  into  three  sections. 
In  §2  we  summarize  the  successful  test  run  on  the  N7.TJ  UNIVAC  of 
a  known  solution  of  equations  (3).  The  results  have  been  stated 
in  a  previous  report  IMM-NYU  2I4.7  (1958)  by  D.  A.  Levine  in  which 
the  numerical  method  of  solution  of  the  Geostrophic  Conservation 
Equation  (1)  is  described.   In  §3  we  describe  the  first  series 
of  2k   hour  runs  (hourly  time  step)  by  the  Total  Flow  Method, 
equations  (1),  using  initial  data  for  Hurricane  lone,  1955.  Al- 
though the  boundary  data  was  held  fixed  in  this  test  series  in 
order  to  simplify  the  machine  coding,  a  partially  successful 
prediction  was  made.   In  iij.  some  partial  runs  by  the  Point  Vor- 
tex Method,  equations  (3),  are  described  and  some  of  the  diffi- 
culties encountered  are  discussed. 


2.   NUMERICAL  SOLUTION  OP  AN  EXACT  SOLUTION 

The  machine  code  described  in  IMM-NYU  2l|7  was  tested 
against  an  exact  solution  of  the  goostrophic  conservation  equa- 
tion (1).  This  test  served  four  main  purposes:   1)   to  test  the 
functioning  of  the  code,  2)   to  verify  estimates  of  error  of 
the  numerical  scheme,   3)   to  establish  the  validity  of  the 
exact  singular  solution  and  l\.)      to  test  the  proper  scaling  of 
physical  quantities. 

The  exact  solution  used  for  the  test  is 

*  =  V*l  =  -  Ttc'V  K.  |r-rQ|  )+Ae-k7  (fca) 

which  satisfies  (1).   \|/,  =  Ae"K7  is  time-independent.  The  mo- 
tion of  the  vortex  by  (3d,e)  is 

-_°  =  „  Ae"  *<>  (Ub) 

dT 

dy 

-£  =  0  (I^c) 

where  we  have  put  the  geostrophic  point  vortex  initially  at  the 
origin.   (Ub,c)  states  that  the  vortex  simply  moves  uniformly 
from  west  to  east; 

xo(t:)  =   A  (I|d) 

y0('l)  =  0  .  ike) 

In  the  numerical  solution,  the  singular  KQ(  lc  |r-rQ  | )  was  trunca- 
ted so  that  the  maximum  tangential  velocity  was  120  mph: 
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The  constants  chosen  were  y/2ti:  =  lCPmiles  hr"  and  K=3»65xlO"^ 
miles"  }  then  the  truncation  or  cut-off  radiu3  r  =  82  miles. 
The  constant  A  was  chosen  so  that  the  translational  velocity  of 
the  point  vortex  KA  =  20  mph.  2h   hourly  time  steps  and  a.  mesh 
width  of  200  miles  was  used  for  a  20  by  20  grid.  The  test  re- 
sults were  remarkably  good: 

1)  The  point  vortex  moved  due  east  at  exactly  20  mph  (to 
six  significant  figures) • 

2)  The  regular  flow  field  \|r,  remained  time-independent  (to 
six  significant  figures). 

The  exact  solution  (J+a)  is  a  special  case  of  a  more  general 
solution  where  the  regular  field  ty^  is 

^  (x,y,  t)  =  f  F(9,'f)exp<+  K.U  cos  e+y  sin  0)  jdQ+constant  ,  (6) 
o  I  J 

with  P(©,^)  an  arbitrary  function.   (6)  may  be  useful  for  addi- 
tional computational  tests. 


3.   NUMERICAL  SOLUTION  BY  THE  TOTAL  PLOW  METHOD 

In  the  total  flow  method  we  do  not  distinguish  between  the 
hurricane  and  the  rest  of  the  flow  field  but  introduce  the  given 
initial  and  boundary  data  directly  into  (1).   It  is  conceivable 
that  if  the  hurricane  is  too  localized  then  it  will  not  be  re- 
solvable by  the  grid  mesh  size  U3ed  (200  nautical  miles).  How- 
ever, we  have  found,  using  initial  data  (500  mb.)  for  Hurricane 
lone,  September,  1955  obtained  from  the  U.  S.  Weather  Bureau, 
that  the  hurricane  does  not  get  lost  in  the  grid  mesh,  although 
the  calculated  closed  contour  lines  around  the  hurricane  do  not 
remain  circular.   In  the  point  vortex  method  described  in  §[)., 
the  vortex  is  necessarily  circular.   Possibly,  with  better  ini- 
tial data  in  the  neighborhood  of  the  hurricane  a  more  nearly 
circular  shape  might  be  retained,  assuming  of  course  that  the 
actual  hurricane  contour  lines  are  approximately  circular. 

Partial  success  was  attained  in  predicting  the  path  of  Hur- 
ricane lone  over  a  2k-   hour  period,  even  though  the  boundary  data 
was  held  fixed  in  order  to  simplify  the  coding.  The  direction 
of  motion  was  simulated  and  the  final  position  (+2lf.  hours)  which 
was  more  than  300  miles  northwest  of  the  initial  position  was  in 
error  by  50-100  miles.  The  non-circularity  of  the  calculated 
contours  does  not  allow  a  better  estimate  of  position.  Although 
the  2l|.  hour  position  was  in  our  opinion  well-predicted  consider- 
ing the  possible  inaccuracies  of  initial  and  boundary  data,  the 
intermediate  hourly  time-steps  show  that  the  calculated  motion 
took  a  large  "jump"  in  the  first  time-step  and  the  motion  there- 
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after  was  extremely  slow.   It  is  not  clear  at  present  whether 
this  behavior  is  due  to  improper  initial  or  boundary  data  or 
both  and  future  tests  are  being  set  up  to  determine  the  causes. 
We  are  more  inclined  to  place  the  blame  on  improper  initial  data. 
Hence  the  following  smoothing  procedure  is  now  being  introduced 
into  the  code: 

1)  Prom  the  initial  data  which  we  call  \Jr,  we  calculate  the 
distribution  F(x,y)  =  (A,-iv  H  a^  each  grid  point. 

2)  Using  this  F(x,y)  we  solve  for  the  smoothed  initial 
data  which  we  call  i|/  by  integration,  (/\,-  K.  H  =  P»  We  expect 
this  smoothing  procedure  to  eliminate  the  abrupt  jump  in  the  nu- 
merical solution.  Also  time-dependent  boundary  data  will  be  used 
henceforth. 

An  indication  of  the  rather  tentative  nature  of  our  test 
runs  to  date  is  illustrated  by  our  uncertainty  concerning  a  good 
estimate  of  the  parameter  K.   Our  tests  so  far  have  indicated 
an  apparent  insensitivity  of  the  numerical  calculations  to 
values  of  K.  differing  by  a  factor  of  two  or  three.  This  insen- 
sitivity seems  very  questionable. 


k*      NUMERICAL  SOLUTION  BY  THE  POINT  VORTEX  METHOD 

In  the  point  vortex  method  vie  carry  out  a  separation  (see 
equation  3a)  of  the  symmetrical  hurricane  field  ty  from  the  re- 
maining regular  field  \|r.  from  the  same  initial  data  (500  mb) 
used  in  §3  for  Hurricane  lone,  September  1955 •  Since  we  have 
not  been  able  to  obtain  reliable  local  velocity  distributions  in 
the  neighborhood  of  the  hurricane  eye,  we  have  encountered  dif- 
ficulty in  getting  a  good  estimate  of  the  hurricane  strength  y 
(see  equation  2a) .   One  method  which  we  have  used  to  obtain  y  is 
by  extrapolation  of  the  distribution  of  the  quantity  r(\|/  )   to 
the  origin, 

-I-  =  Urn  r(y  )  (7) 

2tt   r->0    °  r 

where  r  is  the  radial  distance  from  the  eye.  However,  none  of 
the  available  data  gives  us  much  confidence  in  the  calculated 
values  of  y.  A  general  criterion  which  we  have  used  in  estima- 
ting y  is  to  pick  the  strength  so  that  the  regular  field  \|/,  =  ty-ty 
has  smooth,  nominally-parallel  contours  at  and  near  the  eye  with 
spacing  and  direction  approximating  the  estimated  initial  velo- 
city of  the  hurricane  eye.  Even  if  reliable  data  were  available, 
precaution  must  be  taken  in  separating  \|r  from  the  total  \|/  in 
order  to  avoid  unrealistic  residual  kinks  in  the  contours  of  \|/, 

at  and  near  the  vortex  point  (x  ,y  ) .   In  the  first  test  run 

oo 

which  we  made  by  the  point  vortex  method,  we  apparently  encoun- 
tered such  difficulties.   Unfortunately,  the  coding  for  this  par- 
ticular run  was  set  up  to  yield  the  time-wise  vortex  position 
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10 
(x  ,y  )  only  and  not  the  spacial  \K  distribution  so  that  we  do 
not  have  complete  confirmation  of  our  difficulties.   This  will 
be  remedied  in  future  test  runs.   In  addition,  both  the  smooth- 
ing procedure  for  the  initial  data  (for  \|/,)  and  time-dependent 
boundary  data  will  be  used  in  the  future  (of.  §3  for  total  \|/)  . 
An  outline  of  the  numerical  method  of  solution  by  the  point  vor- 
tex method  is  described  in  IIIM-NYU  2l+7  (1950). 
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